A microfluidic device for investigating crystal nucleation kinetics 
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We have developed an original setup using microfluidic tools allowing one to produce continuously 
monodisperse microreactors (~ 100 nL), and to control their temperatures as they flow in the 
microdevice. With a specific microchannels geometry, we are able to apply large temperature 
quenches to droplets containing a KNO3 solution (up to 50° C in 10 s), and then to follow nucleation 
kinetics at high supersaturations. By measuring the probability of crystal presence in the droplets 
as a function of time, we estimate the nucleation rate for different supersaturations, and confront 
our results to the classical nucleation theory. 



I. INTRODUCTION 

Nucleation kinetics plays a fundamental role in crystal- 
lization processes. In particular, it is responsible for the 
final properties of the crystals: size, size distribution and 
polymorphs. Therefore, the control of nucleation kinetics 
is a real challenge in many scientific fields likeproteomics, 
chemical engineering, 

pharmacology... (HUES,!!- The 

nucleation rate J, the number of nuclei produced by unit 
of time and volume, is the relevant parameter charac- 
terizing nucleation kinetics. In the classical nucleation 
theory, J is given by: 



J = AS cxp - 



B 



{\ogsy 



(i) 



where A is a kinetic factor, S the supersaturation of the 
solution, and B depends on the interfacial tension 7 be- 
tween the nuclei and the supersaturated solution 0,0, SI- 
S' corresponds to the ratio of the activities of the so- 
lution in the supersaturated state and at equilibrium. 
The above expression of J holds for both homogeneous 
and heterogeneous nucleation, only the values of A and 
B differ. The classical nucleation theory accounts well 
for many experimental facts, and some experiments have 
yielded values of A and B which are close to theoreti- 
cal predictions for homogeneous nucleation in the case 
of precipitation S Toll, solidification (polypropylene, n- 
octadecane. . . ) [11. Il2l|. However, homogeneous nucle- 
ation has not yet been observed to our knowledge for 
crystallization of solutes, due to the unavoidable pres- 
ence of impurities. Moreover, it seems nowadays that 
the classical nucleation theory is too much simplified to 
describe the nucleation process in solution. For instance, 
it does not consider the desolvatation energy, and several 
authors propose new mechanisms involving several steps 
during the nucleation O, [3, [H, H] • 

One classical method to estimate the nucleation rate 
J consists of the determination of the induction time U 
which is the time elapsing between the creation of the su- 
persaturation (e.g. obtained after a temperature quench) 



and the first appearance of particles by some detection 
device. One has: 



t% — tn tg 7 



(2) 



where t n is the time required for the nucleation, and t g 
the growth time of the nucleus to the detectable dimen- 
sion. However, it is difficult to determine J from mea- 
surements since it depends on the sensitivity of the de- 
tection device (growth kinetics of the crystals has to be 
known) [§]. Moreover, in this type of experiments, usu- 
ally carried out in macroscopic volumes, impurities play 
a significant role, a rapid temperature quench is difficult 
to achieve, and mixing can have an effect [l7l ]. 

Actually, the most reliable way to measure crystal nu- 
cleation kinetics is the droplets method EE OH H3, 
l2ll . |22| . The principle is to divide a volume of solution 
in a large number of small independant reactors, for ex- 
ample droplets in suspension in an inert oil. In this case, 
one may observe homogeneous nucleation if the number 
of droplets is larger than the impurities initially present 
in the solution. After supersaturation is reached (gen- 
erally by cooling) , one measures the fraction of droplets 
which contain a crystal as a function of time. When 
the droplets volume V is small enough, the nucleation 
time t n is large compared to t g , and so only mononuclear 
nucleation occurs [23|. In that case, nucleation being a 
stochastic process, the probability P that a droplet con- 
tains a crystal evolves as: 



P(t) = 1 - cxp(-cjt) 



(3) 



'Electronic address: philippe.laval-exterieur@eu.rhodia.com 



when the solution is instantaneously supersaturated at 
t = [6j. w = JV corresponds to the nucleation fre- 
quency, namely the number of nuclei formed per unit 
of time. With this method, the simple observation of a 
large number of drops enables one to reach an accurate 
statistics on P and so to determine the nucleation rate. 

Nevertheless, the droplets method raises some exper- 
imental difficulties. Indeed, droplets are generally pro- 
duced by emulsification of the solution in oil using sur- 
factants, and are never perfectly monodisperse. Thus, it 
is necessary to measure their polydispersity to determine 
P(t), and surfactants may induce heterogeneous nucle- 
ation [241. It has also been shown that nucleation events 
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in drops may not be independant in concentrated emul- 
sions [25j. Moreover, detection of the fraction of crys- 
tallized drops is often indirect, and rapid temperature 
quenches are difficult to achieve. 

In this paper, we present a new microfluidic tool suit- 
able for kinetic studies of crystal nucleation. This de- 
vice, based on the droplets method, allows us to over- 
come most of the difficulties described above. Indeed, 
it enables us to produce monodisperse droplets of solu- 
tion without surfactants, to apply very rapid temperature 
quenches, and to measure directly and continuously the 
proportion of crystallized droplets. In a few words, in the 
device, monodisperse droplets of the solution are contin- 
uously produced in an oil phase and convected along a 
microfluidic channel. During the flow, the droplets are 
quickly cooled down (typically 10 s) to a temperature 
below their solubility temperature so that nucleation can 
occur. In the microchannel, the time of flow t of the 
droplets is related to the distance d they cover along the 
channel by t — d/U, with U the droplets velocity. As 
a consequence, we can perform stationary measurement 
of a kinetic process along the flow [2(|. Thus, to deter- 
mine the nucleation rate, we just need to measure the 
probability P at different positions along the channel on 
a large number of droplets to obtain a satisfactory statis- 
tics. In the present work, we use this microfluidic tool 
to measure the nucleation rate of a solution of potassium 
nitrate (KNO3) in water for different supersaturations. 
We show that the device is suitable for rapid measure- 
ments of J as a function of S, and, using the classical 
nucleation theory, to estimate physical properties such 
as the interfacial tension 7. 
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FIG. 1: Design of the microdevice. The oil and aqueous 
phases are injected respectively in inlets 1 and 2. Oil can also 
be injected in inlet 3 to increase the velocity of the droplets. 
The top-right image shows the formation of the droplets. The 
two dashed areas are temperature controlled at T\ and T2. 
The dotted area in the middle of the serpentine corresponds 
to the observation zone through a stereo microscope. 



II. EXPERIMENTAL SECTION 

A. Microfluidic experiments 

Microfluidics refers to systems used to manipulate liq- 
uids in channels with typical dimension of 100 fim (see 
[2H [28| and references therein) . In specific channels ge- 
ometries, a two-phases flow can produce monodisperse 
droplets [2§|. The design of the microdevice is shown 
on Fig. [TJ It is composed of three inlets and one out- 
let. In inlets 1 and 3, viscous silicone oil (Rhodorsil 
500 cSt) is injected at constant flow rate ranging be- 
tween 1 and 4 mL hr _1 . The aqueous solution used 
for the study is previously prepared in a beaker. It was 
made of 83.6 g of KNO3 (Normapur Merck) added to 
100 g of deionized water. The solubility temperature T s 
of this solution is 50°C Q. It is injected in inlet 2 at 
flow rate Q of about 500 /iL hr _1 . The liquids are in- 
jected with syringe pumps (PHD 2000 infusion Harvard 
Apparatus). At the intersection between the oil and the 
aqueous streams, monodisperse droplets are continuously 
produced (see up right insert on Fig. [TJ . It is important 
to note that to avoid nucleation at the surfaces of the 
droplets we choose silicone oil which does not wet KNO3 
crystals. The droplets volume V is controlled by the ra- 
tio of oil (from inlet 1) to aqueous phase flow rates and 
measured by V = Qj /, where / is the droplets produc- 
tion frequency (close to 0.3 droplet per second in our 
device) . With the geometry of the microchannels used in 
these experiments, the droplets volume could be tuned 
from 50 to 200 nL. The monodispersity of the droplets 
was determined by image analyses measuring the time td 
between the production of two successive droplets. Since 
Q = V/td, and Q is constant, variations of td are directly 
related to the variations of V. We find a volume poly- 
dispersity of a few percents. When the flow rates ratio is 
changed to tune the volume of the droplets, the total flow 
rate of liquids can be modified and as a consequence, the 
velocity of the droplets is also changed. To balance these 
variations of the total flow rate, and keep constant the 
velocity of the droplets, silicone oil is injected through 
inlet 3. 



B. Microfabrication 

The device is fabricated in poly(dimethylsiloxane) 
(PDMS) by using soft-lithographic techniques [3p| . 
PDMS (Silicone Elastomer Base, Sylgard 184; Dow Corn- 
ing) was molded on master fabricated on a silicon wafer 
(3-Inch-Si- Wafer; Siegert Consulting e.k.) using a neg- 
ative photoresist (SU-8 2100; MicroChem). Microchan- 
nels of 500 /jm high were used. To make molds of such 
heights, we spin successively two 250 /mi thick SU-8 lay- 
ers on the wafer. After each spincoating process, the 
wafer is soft-baked for 10 min at 65° C and for 60 min 
at 95° C. Photolithography was used to define negative 
images of the microfluidic channels. Finally the wafer 
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was hard-baked for 25 min at 95° C and developed (SU-8 
Developer; MicroChem) . 

In order to have all channel walls made of the same 
material, the PDMS-molded channels were sealed with a 
silicon wafer previously covered by a thin layer of cured 
PDMS of about 50 /j,m. To achieve this sealing, a mixture 
of 80% wt. in PDMS and 20% wt. in curing agent (Cur- 
ing Agent Silicone Elastomer, Sylgard 184; Dow Corning) 
was molded on the SU-8 master described above at 65° C 
for 25 min. At the same time, a mixture of 95% wt. in 
PDMS and 5% wt. in curing agent was spun on the wafer 
to produce a thin layer and was baked at 65° C for 40 min 
|3ll |. The device is then peeled off the mold and holes for 
the inlets and outlets (1/32 in. o.d.) were punched into 
the material. Eventually, the device was placed with the 
wafer covered by the layer of PDMS in a UV ozone appa- 
ratus (UVO Cleaner, Model 144AX; Jelight) and exposed 
to irradiation for 2 min, before sealing at 65° C for at least 
one day. 



C. Observation 

To study the nucleation kinetics in the microdevice, we 
observe all the dotted area on Fig. Q] with a stereo mi- 
croscope (SZX12; Olympus) and a CCD camera (C4742- 
95; Hamamatsu). Thanks to the serpentine shape of the 
channel, the observation of this area allows us to access 
simultaneously different times of the kinetics. A typi- 
cal image of this area is shown on Fig. [2l Such images 
are used to measure the velocity U of the droplets along 
the channel. More precisely, to determine U, we make a 
movie of the observation area during the time required 
for the passage of about 30 droplets through each sec- 
tion. An image processing program, searching the max- 
imum of auto-correlation between two successive images 
of the movie, measures automatically the velocity of each 
droplet. Then, the average of the velocities U is calcu- 
lated for each section of the serpentine, that is to say, for 
different distances along the microchannel. These mea- 
surements, which show on Fig. a slight variation of 
the velocity (less than 10%) probably due to a small in- 
homogeneity of the channel height, enable us to achieve 
an accurate estimation of the time of flow of the droplets. 

The detection of KNO3 crystals formed in the droplets 
is based on their birefringent properties. To measure the 
proportion of droplets that contain a crystal as a func- 
tion of time, we make another movie of the observation 
area but under crossed polarizers and during the time 
t m corresponding to the passage of about 150 droplets. 
Because of their birefringence, we only see crystals (see 
Fig. EJd) . An image processing program counts automati- 
cally the numbers N of crystals that go through each ser- 
pentine section during t m . From the droplets production 
frequency /, we then determine the total number No of 
droplets that go through each section during the movie. 
For an experiment in which N events are observed in N 
trials, the probability to observe an event is P — N/Nq. 
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FIG. 2: (a) Typical image of a movie of the observation area 
in the serpentine. Such images enable us to measure the ve- 
locity of the droplets along the channel, (b) Image of the 
same area obtained under crossed polarizers. Only birefrin- 
gent objects such as crystals of KNO3 can be observed, (c) 
Velocity U of the droplets as a function of their position along 
the microchannel. 



The error on the determination of P is given by the stan- 
dard deviation of the binomial distribution function of P 
in the case of independant events. Experimentally, P is 
determined from a statistics on about 150 droplets, the 
results being the same as for experiments carried out on 
500 droplets. 



D. Temperature control 

To avoid any crystallization before the droplets for- 
mation, the syringe containing the solution and the cor- 
responding tubing are heated at about 60° C with two 
flexible heaters (Minco) and, as shown on Fig. [H the 
inlets area is heated at a temperature T\ also with a 
heater (Minco) placed on the back of the silicon wafer. 
Each heater is controlled with a temperature controller 
(Minco). The serpentine area is kept at temperature Ti 
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using a thermostated support placed underneath. The 
temperature of the support is regulated with a water cir- 
culation and a cryostat (F25; Julabo). 

For the temperature measurements in the device, thin 
thermocouples (type K, 0.5 mm o.d.; Thermocoax) are 
inserted through the PDMS layer down to the silicon 
wafer. These thermocouples enable us to measure the 
temperature close to the microchannel. Figure dis- 
playing a typical temperature profile, shows the local 
heating of the inlets area and the homogeneous temper- 
ature of the serpentine one. We have measured the tem- 
peratures at different positions along the channel, marked 
by the white spots on Fig. for five different minimal 
temperatures Ti of the serpentine area. These temper- 
ature profiles are plotted on Fig. [3Jd. They demonstrate 
that we can apply a large temperature gradient between 
the inlets area and the beginning of the serpentine (typ- 
ically up to 50°C, from 60 to 10°C), maintaining the 
serpentine area at a temperature almost constant, with 
maximal variations of 1 to 2°C. 

The temperature of the liquid flowing along the chan- 
nel has been measured with thermocouples (type K, 
0.5 mm o.d.; Thermocoax) drove into the microchannel. 
We have observed that, anywhere along the channel, the 
liquid is at the same temperature as the wafer. This is 
mainly due to the small velocity of the droplets (typically 
4 mm s _1 ), the small height of the microchannel, the 
forced convection of liquid in the droplets [32], and the 
high thermal conductivity of the silicon wafer. As a con- 
sequence, as soon as the droplets are produced, they are 
quenched down to the serpentine temperature (typically 
50°C in about 10 s) which allows us to access kinetics 
with characteristic times ranging from 10 s to the total 
residence time of about 150 s (for U ~ 4 mm s _1 ). 



III. RESULTS AND DISCUSSION 
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We now discuss the results obtained for the KNO3 
solution whose temperature solubility is 50° C (see ex- 
perimental section). The droplets are formed at 60° C 
and cooled down to different temperatures T2. For 
T2 > 25° C, no nucleation events occur when droplets 
flow along the channel and for T2 < 10°C, crystals sys- 
tematically appear during the temperature quench. Fig- 
ure [4^ presents the temporal evolution of 1 — P for three 
intermediate temperatures T2, where P is the probabil- 
ity that a drop contains a crystal. The origin of the time 
axis is set at entrance of the serpentine. As expected, 
1 — P decreases with time and the kinetics of crystal nu- 
cleation becomes faster as T2 decreases (Eq. ([3])). Note 
that, at the investigated range of temperatures, growth 
kinetics (experimentally measured of the order of 1 s) is 
very fast compared to the nucleation kinetics (typically 
50 s) so we can assume that the time at which a crystal 
is first observed corresponds to the nucleation event. As 
shown by Eq. , 1 — P is expected to follow an exponen- 
tial decay characterized by the nucleation frequency u>. 



FIG. 3: (a) Temperature profile of the device. The gray scale 
is linear with the temperature. The hot area is in white and 
the cold one is in dark grey. The hatched disk corresponds 
to the silicon wafer, (b) Temperature of the droplets flowing 
along the channel for different temperatures T2 of the cooled 
area. The positions along the channel where the temperatures 
are measured, are marked by white spots on image (a). 



However, the results displayed in Fig. HJd, can not be ac- 
counted by simple exponential behaviours but rather by 
1 — P = cxp(— Lut+e). Such offsets at t = s have already 
been reported and can be explained by t he p resence of im- 
purities in some of the droplets [13, EH, [33| . Indeed, they 
may provoke heterogeneous nucleation with faster kinet- 
ics. Therefore we observe a fraction of droplets which 
contain crystals before the first observed channel (« 25 s 
in our device). Surprisingly, slightly positive offsets, cor- 
responding to delay times for the kinetics, have been sys- 
tematically observed at temperatures higher than 20° C, 
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FIG. 4: (a) 1 — P and the droplets temperature (dotted plots) 
as a function of time. P is the probability to have a crystal 
in a droplet, (b) log(l — P) as a function of time. These 
measurements have been obtained with droplets of 100 nL and 
for three different temperatures of the cooled area (bullets: 
13° C, squares: 16° C and triangles: 20° C). The time origin 
corresponds to the entrance of the droplets in the serpentine. 



temperature and therefore, that nucleation rates can be 
estimated by J = w/V. Even if the range of volume 
variation is not very large, these results indicate that lj 
is proportional to the number of nucleation sites inside 
the droplets and that nucleation occurs in the droplet 
volume. 




where probabilities P do not exceed 5% at t = 150 S. This 
seems to indicate that there is no fast kinetics induced 
by impurities as above. In this case, the observed delay 
times may be due to slight variations of the temperature 
along the serpentine which can play an important role at 
low probabilities. 

Despite the observed offsets, the nucleation frequen- 
cies uj are estimated by fitting the data according to 
1 — P = cxp(— u)t + e). lj are almost constant over the 
investigated range of time. Figure [5^ presents the nucle- 
ation frequency w as a function of the droplets volume V 
for different temperatures T 2 (in all these experiments, 
the velocity of the droplets is almost the same). Our 
results suggest that co is proportional to V at a given 



FIG. 5: (a) Nucleation frequency u as a function of droplet 
volume V for different temperatures of the cooled area (bul- 
lets: 13°C, squares: 16°C, diamonds: 18°C and triangles: 
20° C). The solid lines are the linear fits of the plots for the 
different temperatures, (b) Probability p to have n droplets 
without crystal between two droplets containing a crystal, for 
P = 0.41. The bars correspond to the experimental data 
obtained for 200 droplets and the solid line is the geometrical 
lawp(n) = P{1-P)". 



In the droplets method, it is important to verify that 
the presence of a crystal in a droplet does not influence 
the nucleation in others (25| (e.g. contamination of the 
channel in our device). When the nucleation events are 
independent, the probability p to observe n droplets with- 
out crystal between two droplets containing a crystal, 



must follow the geometrical law 



and in the case of spherical nuclei, B is given by 



p(n) = P(l - py 



(4) 
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We have measured the probability p(n) for different prob- 
abilities P in several experiments. On Fig. [5b, we 
compare the experimental results obtained at 16° C for 
P = 0.41 with the corresponding geometrical law. Good 
agreements between the two plots are systematically ob- 
served, especially when the statistics is carried out on a 
large number of droplets. We conclude that the presence 
of a crystal in a droplet does not induce nucleation in 
others, and our device is therefore well-suited to follow 
independent nucleation events. 

By changing the temperature T%, we could measure 
the nucleation rate J for different supersaturations S. 
To calculate S, we assumed that S w (C„/C e(? ) 2 , where 
Cin and C eq are respectively the initial and equilibrium 
concentrations of KNO3 (despite numerous literature re- 
searches, no experimental data could be found concerning 
the ions activity at the investigated metastable states). 

In the performed experiments, S has been varied in 
the range of to 13. Below S = 5, no crystals nucleate 
within the investigated time (150 s) and above S = 13, 
crystals systematically appear during the temperature 
quench (10 s). Figure [6] presents log J/S as a function 
of l/(log5) 2 . The error bars take into account the small 
temperature variations along the serpentine. Within the 
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FIG. 6: (Plot of log J/S as a function of l/(log S) 2 . J is the 
nucleation rate (m -3 s _1 ), S the supersaturation. The solid 
line is the linear regression of the data. 



experimental uncertainties, Eq. |TJ well fits our data with 
A ?s 3.10 7 m -3 s _1 and B w 12. In this case, we have 
supposed that both A and B do not depend on temper- 
ature. In accordance with the classical nucleation theory 



167T7 v 

3(fcT) 3 



(5) 



where 7 is the interfacial tension between the nucleus and 
the solution, v the molecular volume in the crystal, and k 
the Boltzman constant. Supposing that 7 is proportional 
to T (i.e. B does not depend on T), we found 7/T = 
6.6 10~ 2 mN m" 1 K" 1 with v = 8 10" 29 m" 3 0. In the 
investigated range of T, 7 ps 19 ±1 mN m . Note that if 
we suppose 7 does not depend on T, Eq. (TJ and Eq. (j^J) 
give also a correct fit of the data with 7 pb 20 mN m . 

In the case of homogeneous nucleation, theoretical es- 
timations of A are of the order of 10 26 -10 30 m~ 3 s _1 . 
The A we measured (A « 3 10 7 m~ 3 s _1 ) suggests we 
observe heterogeneous nucleation. Moreover, over sev- 
eral series of experiments (corresponding to observation 
of about 6000-8000 droplets), we found a small disper- 
sion of the 7 values (16-26 mN m _1 ) and values of A 
ranging from 3 10 7 to 10 10 m -3 s _1 also suggesting a 
heterogeneous mechanism. Furthermore, since the typ- 
ical concentration of impurities is about 10 -10 s mL _1 
according to 0, each droplet of 100 nL contains about 
10 2 to 10 4 impurities susceptible for acting for heteroge- 
neous nucleation. 

Nevertheless, it is difficult to tell if we observed ho- 
mogeneous or heterogeneous nucleation just from com- 
parisons with the classical nucleation theory. Indeed, up 
to now, no experimental values of A as large as 10 26 - 
10 30 m~ 3 s _1 have been measured for the crystallization 
of solutes (except for the case of precipitation where high 
S can be reached), and the classical nucleation theory 
may not take into account the entire complexity of solute 
crystallization 0, 0, H, 0, EH- For exemple, according 
to the classical nucleation theory, the kinetic factor A is 
proportional to the concentration of nucleation sites in 
the solution Q. As a consequence, we can write: 



A 



HO 



.4 



HE 



Co 



(6) 



where Ano and ^4he are the kinetic factors in the case of 
homogeneous and heterogeneous nucleation respectively, 
Co is the molecule concentration in the solution, C, is the 
concentration of heterogeneous nucleation sites and a is 
a factor of the order of 1 [1] . With Aho 
Co = 10 28 m~ 3 and A^e = 10 7 m -3 s~" m our case, 
we have C ~ 10 9 m -3 . Since the droplets volume is 
100 nL, this results would suggest that there are less than 
one heterogeneous nucleation site per droplet, which is 
contradictory with our previous conclusion. 

To bring new experimental insights, one should in- 
crease the supersaturation 0, 0| in the solutes crystal- 
lization experiments. In our device, all crystals nucleate 
during the temperature quench at S > 13 for droplets 
volume of 100 nL. To reach higher supersaturations and 
measure the corresponding kinetics, we should decrease 
significantly the droplets volume down to a few nL, such 
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tiny droplets would allow us to access homogeneous nu- 
cleation [fj]. 

IV. CONCLUSION 

To conclude, we have performed original microflu- 
idic experiments suitable for nucleation kinetics measure- 
ments. Indeed, our microdevice enables us to overcome 
the difficulties encountered with the classical droplets 
method. More precisely, we have managed with our de- 
vice to apply rapid and large temperature quenches to 
monodisperse microreactors, to transport solid particles 
in droplets, and to monitor continuously with direct ob- 
servations, the fraction of crystallized droplets as a func- 
tion of time. Such experiments yield rapid statistical 
measurements of the nucleation kinetics and therefore 



give access to value of nucleation rates with a low liquid 
consumption (200 droplets w 20 fxh). 

We also believe our system is well-adapted to follow 
kinetics of any temperature dependent processes. With 
only a few enhancements of the experimental device 
(automation of several liquids injection to vary solvent 
and solutes concentration) , high-throughput screening of 
crystallization conditions could be carried out. 
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